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Abstract 

The searches for H — > WW events at the LHC use data driven techniques for estimating 
the qq — > WW background, by normalizing the background cross section to data in a 
control region. We investigate the possibility that new physics sources which mainly 
contribute to the control region lead to an overestimate of Standard Model backgrounds 
to the Higgs boson signal and, thus, to an underestimate of the H — > WW signal. A 
supersymmetric scenario with heavy squarks and gluinos but charginos in the 200 to 
300 GeV region and somewhat lighter sleptons can lead to such a situation. 
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1 Introduction 



Among the Higgs search channels at the LHC, the decay into weak boson pairs is particularly 
important, because the tree level couplings of the Higgs boson to WW and ZZ allow to 
positively identify a scalar state as being associated with the vacuum expectation value 
which breaks the electroweak gauge symmetry. A Higgs boson produced via gluon fusion 
with subsequent decay into leptonically decaying W bosons is especially searched for, at the 
present time, because this channel is expected to produce the largest sample of H — > VV 
events at the LHC. However, for a Higgs boson mass around 125 GeV, it does suffer from 
substantial backgrounds, the dominant one being qq — > W + W~ , followed by leptonic decay 
of the W bosons. In the present experimental analyses, this contribution to the background is 
estimated with data driven methods [1-3]. A control region of the phase space, namely events 
at large dilepton invariant mass m«, is defined where no signal events are to be expected. 
The size of the background is measured in this control region and the corresponding value 
in the signal region, at small m«, is then extrapolated using the shape of the distribution 
determined by a Monte Carlo simulation for Standard Model (SM) W pair production. 

In this paper we investigate whether new physics contributions can seriously compromise 
this data driven background determination. Specifically, can new physics events, which typ- 
ically are hard and thus tend to preferentially populate the large mm region, significantly 
enhance the size of the background measured in the control region, while contributing rel- 
atively less than the SM W pair production to the signal region? This would lead to an 
over-estimate of the background in the signal region and thus to a significant underestimate 
of the Higgs signal. Since the presently measured H — > WW rates are indeed below SM 
expectations, for uih ~ 125 GeV [2,3], such scenario deserves serious investigation. 

We will focus on processes arising within the Minimal Supersymmetric Standard Model 
(MSSM) as an example for processes induced by physics beyond the Standard Model (BSM 
physics). SUSY processes contributing as background have been discussed before with re- 
spect to calibration processes [4] and to Higgs boson searches within vector boson fusion 
where the Higgs boson decays into two tau leptons or two W bosons [5]. Also, in [6] SUSY 
processes as a possible background in the signal region of Higgs boson searches have been 
considered. 

2 Scenario and analysis setup 

Within the MSSM, processes that can contribute significantly to the signature of W boson 
pair production involve the production of charginos, sleptons and neutralinos and therefore 
the parameters that influence their masses are the most relevant ones. 

Our scenario is based on the light slepton scenario described in [5], which is tuned for 
sizable chargino and slepton pair production. The tau slepton masses are chosen to be 
above the light chargino mass while the sleptons of the first two generations are lighter 
than the light chargino. As a consequence, the main decay channels of the charginos are 
via a slepton of the first two generations and a corresponding lepton, where the slepton 
decays directly into a lepton and the lightest neutralino (which is assumed to be the lightest 
supersymmetric particle (LSP)). We increase the chargino mass by modifying the soft SUSY 
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breaking parameters M 2 and to be in agreement with the ATLAS 2 fb~ trilepton 
search [7]. However our scenario shows strong tensions with the CMS trilepton analysis 1 [8] 
which uses the full 5 fb _1 dataset of 2011, if the predicted neutralino properties are taken 
literally. Since production processes with next-to-lightest neutralinos are not important for 
our analysis (see below), we ignore this tension having in mind alternative models that avoid 
the visible leptonic signature of the xl- 2 The masses of squarks of the first two generations 
and of the gluino are not important for the following study and are set to high values. The 
parameters in the Higgs boson sector of the MSSM are assumed to have values which lead to a 
Standard Model type Higgs boson with mass of 124.7 GeV, which is within the experimentally 
allowed Higgs boson mass range and in the mass area where some experimental hints of a 
Higgs boson have been observed [9,10]. Starting from this "base" scenario we vary the 
parameters that control the chargino, slepton and lightest neutralino masses to study their 
influence on the production and decay of the SUSY particles and their contribution to the W 
pair production signature. A scenario with 25% higher LSP and 40% higher slepton masses 
is discussed in more detail. The full parameter settings of the SUSY scenarios can be found 
in Appendix A. The following SUSY production processes give the dominant contributions 
to the signature of the W boson pair production: 



Chargino pair production gives the by far largest contribution, followed by the slepton pair 
production processes. The production of a next-to-lightest neutralino and a chargino gives 
only a very small contribution to the W boson pair production signature, as it produces 
mostly three leptons and the third lepton can be tagged quite well. Therefore this channel 
is not crucial for our results and general BSM scenarios without such a trilepton source 
can evade detection and still show the same behaviour concerning the WW background 
estimation. 

It should be emphasized that the important feature of the scenarios is the existence of 
new physics particles which decay into electrons, muons and invisible particles, which is the 
same particle content in the final state as the one of the W pair production with subsequent 
leptonic decay. In that respect we consider our MSSM scenarios as an example for BSM 
physics. 

The event generation for our analysis is done using Herwig++ 2.5.2 [11] at parton level, 
including parton shower. Within Herwig++, the W pair production is generated at next- 
to-leading order QCD within the POWHEG framework [12]. The SUSY particle pairs are 
generated at leading order [13]. Their cross section is multiplied with an appropriate K- 
factor (K = 1.2), which we obtain from Prospino2 [14]. Spin correlations within production 
and decay in Herwig++ are included as described in [15]. A comparison with a combination 
of MadGraph 5.1.3 [16] and Pythia 6.4 [17] led to comparable results. 

1 A possibility to weaken this tension is to allow a decay of the light chargino and the next-to-lightest 
neutralino into all three lepton flavors in equal parts. However this reduces the effect on the correction 
factor C by roughly 35-40%. 

2 This can be achieved for example with a scenario (typically beyond the MSSM) where the next-to- 
lightest neutralino is much heavier than the light chargino. 



qq -> xT Xi 

Q q -> Xi X°2 



t + r+f T 



(i) 

(2) 
(3) 
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Both ATLAS [2] and CMS [3] have presented a study of the H -»■ WW 2£2v channel 
with the full data set of 2011. We largely use the cuts and methods from the ATLAS analysis, 
because they show distributions in the transverse mass of the W boson pair [18] 

m T = ^ {E e T e + E™ iss ) 2 - |p£ + p™ ss | 2 , with E% = ^\p%\ 2 + m 2 u , EJp iss = \p% iss \ (4) 

in the signal and control regions up to quite high values in Wt- This gives an opportunity 
to check if the effects of BSM physics could be identified in the experiment as an excess of 
particularly hard events. CMS performs a similar analysis, but they only show distributions 
of the invariant lepton pair mass mu and the azimuthal angle A0« between the leptons. 
Those distributions have turned out to be less illuminating for our study. However, as the 
CMS cut selection is similar to the ATLAS one, our results should hold qualitatively for 
CMS as well. 

Our analysis is carried out for the LHC operating at a center-of-mass energy of 7 TeV. 
We use the CT10 parton distribution functions [19] for the POWHEG event samples and 
cteq611 pdfs for the leading order SUSY calculations. As renormalization and factorization 
scales we use the invariant mass of the W boson pair or the SUSY particle pair. For the 
basic event selection we require two oppositely charged leptons £ (electrons or muons), where 
the harder lepton with respect to transverse momentum is labeled £\, the softer one t^. 
The following cuts are applied, largely taken from [2]: 

p TA > 25GeV Vtm > 15GeV 
m eeM > 12GeV m efl > lOGeV (5) 
\m eeM -M z \ > 15GeV \q t \ < 2.5. 

The events are categorized according to the number of visible jets. Jets are clustered using 
the anti-&t algorithm [20] with distance parameter R = 0.4 and the following requirements 
on rapidity rjj and pr t f 

\rjj\ < 4.5 p Tij >25GeV. (6) 

Leptons that are within the R-separation AR < 0.3 of a jet are counted as part of the jet. 
For QCD background suppression in the H — > WW analysis the LHC experiments use the 
quantity E™™^ = E™ lss ■ sinmin(A0, |), where E™ lss is the missing transverse energy of the 
event and A0 is the azimuthal angle between the E™ tss vector and the nearest lepton or jet 
with pt > 25 GeV. The requirement is 

E%$% > 45 GeV for ££ = ee (jifi) E™% > 25 GeV for ££ = e/i . (7) 

The spin-0 nature of the Higgs boson is exploited by demanding [21] 

N>u < 1-8 • (8) 

The signal region in the 0-jet channel is furthermore restricted in the dilepton transverse 
momentum p^ 

p ee( MAt)>45GeV ^> 30 QeV (9) 

and [22] 

mu < 50 GeV . (10) 
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Figure 1: Invariant lepton pair mass mu (left) and transverse mass ttit (right) distributions 
for the "base" scenario of Section 2. The my distribution is calculated including all control 
region cuts of Eqs. (5) - (7), (9) and (12). For the m u plot the m^-cut of Eq. (12) is omitted. 
Both plots show the qq — > WW distribution, the SUSY contributions and their sum. The 
m« plot also shows the qq — > WW result, rescaled by (&$ w + &c ) / a cf W ■> extracted from 
the control region. 

For the 1-jet channel ATLAS uses cuts on the vectorial sum of the Pt of jets, leptons and 
p™ ss , and on the rr invariant mass m TT , calculated in the collinear approximation [23] 

|Pt + Pt + Pt + Pr*"l < 30 GeV \ m rr ~ M z \ > 25 GeV . (11) 

Events with identified b-jets (80% efficiency, 6% mistag) are rejected. The WW control 
regions for the 0-jet and 1-jet bin are defined by omitting the A0« and mee, ma x cuts of 
Eqs. (8) and (10) and requiring a minimal invariant lepton pair mass of 

m u > 80 GeV . (12) 

We do not consider the 2-jet channel, as there is not enough statistics at the moment for 
any conclusions in this channel. Detector effects, efficiencies and hadronization effects have 
been neglected. 

3 Results for SUSY example 

The effects of the SUSY processes of Eqs. (1) - (3) on the full m u range (signal 3 and control 
region) in the 0-jet channel can be seen in Figure 1. The situation for the 1-jet channel is very 
similar, 4 but the BSM effects are much less constrained compared to the 0-jet channel due 
to smaller event numbers. Therefore we focus on the 0-jet bin. Chargino pair production 
accounts for the largest part of the SUSY signal, especially in the control region, while 
slepton pair production has larger effects in the low m« region due to the small assumed 



3 The full signal region cuts also include a cut on Acf>u, Eq. (8), but the effect on the WW and SUSY 
processes is marginal in the m« < 50 GeV region. 

4 Additional jets for the SUSY processes are simulated by the parton shower of Herwig++. 
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Figure 2: Event numbers of the SUSY contributions in the signal region, in the control 
region and in the control region with > 350 GeV and > 440 GeV for varying slepton 
masses of the first two generations. The LSP mass is m x o = 99 GeV in the left plot and 
m x o = 124 GeV in the right plot. The discussed "base" and "worst case" scenarios are 
marked. 



slepton masses. As stated before, the production of a light chargino and the next-to-lightest 
neutralino has only a very small contribution and is therefore not important for our results. 
The relative contribution of the SUSY processes to the signal region is clearly much smaller 
than the contribution to the control region. Therefore this scenario is potentially dangerous 
for the data-driven estimation of the qq — > WW background: If the BSM physics could not 
be identified, the WW prediction in the signal region would be rescaled by 

{a c +a c )/a c , (13) 

where Gq W and aQ USY are the WW and BSM contributions in the control region. In our 
example, this leads to a WW prediction for the signal region which is clearly too high (see 
Figure 1). Furthermore the effect on the shape of the m« and ttit distributions in the signal 
region is too small for a detection of the SUSY contamination. 

In contrast, a closer look at the transverse mass distribution can reveal the BSM physics 
effects of this scenario. As m-r is bounded from below by m« and additional missing trans- 
verse energy results in even larger values of m^, BSM effects with large m« naturally lead 
to contributions at high itit values. This is especially the case in theories with additional 
sources of missing energy like the MSSM. Figure 1 shows an enhancement due to the SUSY 
contributions of more than 100% for my values exceeding 350-400 GeV. ATLAS measured 
41 events with > 350 GeV, with a total background expectation of 48 events, includ- 
ing a WW contribution of 31 events [2]. Therefore a factor of two increase in the U WW 
contribution" is already ruled out with current data. 

Most of the SUSY contributions arise from chargino pair production. As the masses of 
the chargino decay products play an important role for the kinematics of the final state lep- 
tons and for the amount of missing transverse momentum, we vary the soft SUSY breaking 
parameters Mi, M e L and M^l, which govern the LSP and left-handed slepton masses. The 
slepton mass variation also directly modifies the slepton pair production contributions to sig- 
nal and control region. For the identification of a potentially dangerous scenario (concerning 
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the normalization of the WW background with the help of a control region) the following 
constraints have to be fulfilled: 

• The contribution to the signal region has to be as low as possible. 

• The contribution to the control region has to be as large as possible, but still small 
enough to hide in the shape uncertainties of the control region. 

• The part of the control region with high my is strongly constrained by the current 
ATLAS data. Therefore the BSM effect in this region has to be small. 

For the comparison with the ATLAS data from [2], we convert our cross sections into 
expected number of events by normalizing our WW prediction for 4.7 fb -1 with the expected 
number of events from the ATLAS H — > WW study within the 0-j control region. From this 
rescaling we estimate an overall efficiency of 59% for the evolution of showered parton level 
events to reconstructed jets and leptons in the analysis. For the rescaling we also take the 
gg —> WW contribution into account using gg2WW [24] , which is included in the pp — > WW 
background of the ATLAS study. This part is formally of next-to-next-to leading order 
in QCD with respect to the qq — > WW contribution and adds a few percent to the cross 
section [1]. For the rest of our study, the gg — > WW contribution has been neglected. 

With this prescription the number of events shown in Figure 2 is calculated for a LSP 
mass of 99 GeV (left diagram) and 124 GeV (right diagram). For each plot the slepton mass 
is varied up to the chargino mass, bounded from below by the requirement that the lightest 
neutralino has to be the LSP. These event numbers have to be compared with the following 
values for the qq — > WW prediction (our Monte Carlo prediction, scaled with the overall 
efficiency factor 0.59): 

N WW _ one N ww — 454 i\rww _ 9o N ww _n / 14 x 

1 V signal ~ OOD ^ control ~ JV m T >350GeV — ZZ JV m T >440GeV — ' J 

or A^^ 350GeV = 31 and A^^ 440GeV = 11 as taken from the ATLAS trt distribution [2]. 
The discrepancy between those numbers can be ascribed to higher efficiencies for high itlt 
events. However, for our study we can largely eliminate this uncertainty by comparing ratios 
Nmcj Y ' and the accuracy is sufficient for an approximate comparison with ATLAS 
data. 

A slepton mass roughly in the middle between the LSP and chargino masses gives the 
largest contribution to the control region. At the same time, the tail of the tut distribution 
is significantly smaller than for lighter sleptons. Furthermore the larger slepton mass shifts 
the slepton pair production contribution from the signal region to the control region. The 
increased LSP mass in the right plot of Figure 2 leads to less available kinetic energy for 
the decay products and therefore to smaller mu, mx and p™ lss . This further reduces the 
contributions in the high mr bin of the control region. We also considered a larger chargino 
mass as input. This reduces the chargino pair production cross section and therefore leads 
to a reduction of the overall effect, as expected. 

Taking the criteria as described above, we find a scenario which is compatible with the 
ATLAS H — > WW analysis of 2011 data. We refer to it as the "worst case" scenario. As can 
be seen in Figures 2 and 3, this scenario gives very small contributions to the high part 
of the control region, small enough so that they cannot be identified at the moment. At the 
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Figure 3: Invariant lepton pair mass m« (left) and transverse mass ttit (right) distributions 
for the "worst case" scenario of Section 3. The distribution is calculated including all 
control region cuts of Eqs. (5) - (7), (9) and (12). For the m u plot the m u -cut of Eq. (12) 
is omitted. Both plots show the qq — > WW distribution, the SUSY contributions and their 
sum. The mu plot also shows the qq — > WW result, rescaled by {(Jc W + &c USY ) I a 'c W ^ 
extracted from the control region. 



same time the contributions to the signal region are very small and therefore not noticeable, 
although they are at partially higher values of wit than the WW background. 

We now want to quantify the effect of this BSM scenario on the WW background pre- 
diction calculating the factor C by which the expected number of WW events in the signal 
region obtained from the normalization would have to be corrected. The extrapolation for 
the number of events from control to signal region is done using [1] 

jyWW 

N s =^ L N c = a-N c . (15) 

Taking both Standard Model WW production and the BSM effects in the control region into 
account, this leads to a predicted number of background events in the signal region given by 

N norm = Q . (jyWW + pjSUSYj ? ^ 

while the actual contribution is 

N true = N WW + N SUSY 

Therefore the predicted number of events would have to be reduced by 



C 



AT true „WW , „SUSY „WW 



Ajnorm ^WW , ^SUSY ^WW 

i\ s o c -|- a c a q 



where we have replaced ratios of event numbers by our theoretical cross sections. For this 
specific scenario we get 

C = 0.924. (19) 
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For the ATLAS data which enters their fitting procedure and the resulting exclusion 
limits a selection cut of 

0.75 • ran < rriT < tuh (20) 

on the transverse mass is applied. Within this range the relative BSM contribution of our 
"worst case" scenario is even smaller, leading to a larger correction for the extraction of the 
WW background. In this case the number of WW events would have to be reduced by a 
factor of 

C = 0.897. (21) 

Since a SM Higgs signal is about 20% of the overall background, an overestimate of 10% in 
the (dominant) WW background would lead to a very large underestimate in the size of the 
extracted Higgs signal. 



4 Conclusions 

Data driven methods for background determination are extremely useful for reducing the- 
ory errors inherent in QCD predictions for LHC cross sections. However, they do rely on 
assumptions, namely that, apart from the searched signal events, there are no other BSM 
contributions which could affect the search. In this paper we have studied the impact of 
new physics contributions on the estimate of the SM background to the H — > WW — > Ivlv 
search. In the case of W pair production as background to Higgs boson searches, the number 
of events is measured in a high m u control region where no signal events are to be expected. 
Via an extrapolation using Monte Carlo predictions for the SM-shape of the distributions of 
W pair production, the estimate of the background in the softer signal region is determined. 
In general, new physics at high energy scales can enhance the number of events in hard 
control regions while contributing little to a substantially softer signal region. The larger 
measured rate in the control region then can lead to an overestimate of the background in 
the signal region and dilute a potential signal. 

Two example scenarios in the context of the MSSM have been discussed in some detail. 
In the first one, the new physics contributions not only lead to an enhancement of the event 
rates in the control region but also change the shape of distributions (the ttit distribution in 
the case considered) sufficiently to make the extra BSM contributions noticeable. The new 
physics contributions in the second example are more difficult to distinguish from the pure 
SM case and could indeed have been missed in the H — > WW analyses. With enough data, 
one could, of course, see deviations in the hard event distributions, like in the tail of ttit 
distributions for the scenarios at hand. 

The MSSM scenario described above is just one example of BSM physics which might af- 
fect the Higgs search and, as importantly, the measurement of Higgs couplings from measured 
Higgs rates. Such potential BSM contamination should be kept in mind when interpreting 
the Higgs search data within BSM scenarios: the errors on Higgs boson couplings may be 
larger than in a pure SM analysis. 
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A Parameters 

The "base" scenario is determined by the soft SUSY breaking parameters 





M x = 


103.1 GeV 


M eL = M flL = 


134.4 GeV 




M 2 = 


270.1 GeV 


M eR = M llR = 


135.8 GeV 




M 3 = 


1703.7 GeV 


m tL = 


393.6 GeV 




A t = 


-2194.8 GeV 


M rR = 


333.4 GeV 




A b = 


-1907.2 GeV 


M qiL = M q2L = 


1579.8 GeV 




A r = 


-249.4 GeV 


M uR = M cR = 


1524.3 GeV 


A u 


= A C = 


-655.5 GeV 


M dR = M sR = 


1517.7 GeV 


A d 


= A S = 


-821.8 GeV 


M g3 L = 


1201.4 GeV 


A e 


= A» = 


-251.1 GeV 


M tR = 


1019.4 GeV 




M 2 Hd = 


32609 GeV 2 


M bR = 


1257.2 GeV 




M 2 H = 


-169877 GeV 2 


tan/3(M z ) = 


10.0 



(22) 



at the scale Q — 1 TeV and the following Standard Model parameters, with the top mass 
from [25]: 

a~^{M z ) = 127.934 G F = 1.16639 ■ 10" 5 GeV" 2 

a s (M z ) = 0.1172 M z = 91.187 GeV (23) 

M b (M b ) = 4.25 GeV M t = 173.2 GeV . 

These parameters are fed into SUSYHIT [26] for the calculation of the SUSY particle masses 
and branching ratios. The SLHA output file is then used as input for FeynHiggs 2.8.6 [27] 
in order to get precise Higgs boson mass values. The resulting scenario exhibits the following 
features: 

• The squark mass values of all three generations and the gluino mass values (rriq m 
1581 GeV, m Ti = 934 GeV, m~ hi = 1232 GeV, m~ g = 1725 GeV) are beyond current 
exclusion limits [28-31]. 

• The trilinear coupling A t is adjusted according to the maximal mixing scenario [32], 
which yields a Higgs boson mass with values in the vicinity of the experimental hints 
of a Higgs boson [9, 10]. 

• The wino mass parameter M 2 and the Higgs mass parameter m# u are chosen to give 
chargino masses outside the exclusion limits of the ATLAS trilepton search [7]. The 
discussion of the 5 ftT 1 CMS trilepton analysis [8] can be found in Section 2. 
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"base" scenario 


"worst case" scenario 




Mi = 103 GeV 


M e ^ L>R =lA-M^ L>R , 
Mi = 129 GeV 


A.1 

m xi 


98.9 GeV 
260.0 GeV 
260.3 GeV 


124.1 GeV 
260.3 GeV 
260.7 GeV 


™e L = 

m e R = m Ji R 


141.8 GeV 
142.6 GeV 


193.5 GeV 
195.0 GeV 


BR( X +^l+Z £ ) 
BR(xt^W + X °i) 


58.0 % 
41.0 % 
1.0 % 


60.4 % 
37.7 % 
1.9 % 


Bft W2 ^ l L,R l ) 

BR{ x ° 2 ^viis e ) 
BR(xl^x\Z) 
BR( X ° 2 ^x1h ) 


A7 7 % 
^± 1 . 1 /o 

51.4 % 
0.2 % 
0.7 % 


^O.l /o 

53.2 % 
0.3 % 
1.4% 


BR(£i >R ^xW 


100.0% 


100.0% 



Table 1: Masses and branching ratios of interest for our "base" scenario and for a scenario 

with modified Mi and M e ^L tR . Here, I means the combined (s)electron and (s)muon 
channel. For completeness, we also give the mass and the branching ratios of the x\ though 
they are not important for our results. 



• The stau lepton masses are larger than the light chargino mass (m^ = 334 GeV), which 
is a specific feature of the considered scenarios. 

• The mass parameters for the left-handed sleptons of the first two generations M e £ 
and M^l are chosen such that the chargino decay into selectrons and smuons is the 
dominant decay mode. 

The soft SUSY breaking parameters of the "worst case" scenario discussed in Section 3 are 
the same as in the "base" scenario, except for 

Mi = 128.9 GeV M eR = M„ R = 190.1 GeV 

M eL = M flL = 188.2 GeV. 1 ' 

The relevant masses and branching ratios of both scenarios are listed in Table 1. 
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